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ABSTRACT 

We propose a two parameter generalization for the dark energy equation of state 
(EOS) for thawing dark energy models, which includes PNGB, CPL and algerbraic 
thawing models as limiting cases; and confront our model with latest Supernova Type la 
(SNe la) Data from Union 2.1 compilation, latest Observational Hubble Data (OHD), 
Cosmic Microwave Background (CMB) Data from 7 year WMAP results and latest 
Baryon Acoustic Oscillation Data from SDSS Data Release 9 to constrain our parameter 
space. The analysis reveals that thawing dark energy EOS is not unique from the 
observational point of view so far as the standard parameters (dark energy EOS, present 
value of matter density parameter and Hubble parameter) are concerned. In other 
words, different thawing dark energy models are not distinguishable from each other 
with the help of best-fit values of matter density parameter at present epoch, linear 
growth of matter perturbation and the average deceleration parameter. But tuning the 
model parameters does leave its impact on the variation of dark energy EOS as well 
as the model-independent parameters like the statefinder pair and Om3 parameters. 
We are thus led to the conclusion that unlike the standard parameters, the model- 
independent parameters and the variation of EOS serve as model discriminators for 
different thawing dark energy models. 



1. Introduction 

During the past two decades strong observational evidences have come out which have made 
late time acceleration of the universe almost a de facto phenomenon. These evidences have 
been brought forth d la independent astrophysical obs ervations like Supernoy a e Type la (SNIa 



luminosity distance modulus as a function of redshift flA manullah et al. 2010; 



Kessler et al. 20091 : iKowalski et al. 20081 : IPerlmutter et al. 19991 : iRiess et al. 1 998: 



Davis et al. 2007 



Riess et al. 2007 



1 Astroparticle Physics and Cosmology Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 
700064, India; debabrata.adak@saha.ac.in 

2 Astroparticle Physics and Cosmology Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 
700064, India; debasish.majumdar@saha.ac. in 

3 Physics and Applied Mathematics Unit, Indian Statistical Institute, 203 B.T.Road, Kolkata 700108, India; 
supratik@isical. ac.in 



-2- 



Riess et al. 20091: ISuzuki et al. 2003: IWood-Vasev et al. 200711. Observational Hub ble Data (OHD) 
(jGaztanaga et al. 2009 ; I Jimenez et al, 20 02; Sim on et al. 20051: IZhang et al. 201 2j). Cosmic Microwav e 



Background (CM B) Shift Parameter (jKomatsu et al. 2009 ; iKomatsu et al. 2011 



Podariu et al. 2001 



Ratra et al. 19991) an d Baryon Acoustic Oscillation (BAO) Data (jPercival et al. 2010l ; ISanchez et al. 20121 ; 
Samushia et al. 2 009). Attempts have been made to explain accelerated expansion on the ba- 
sis of existence of some exotic fluid, namely dark energy, in huge abundances in the universe . 
Though there exists a lot of dark energy models (see for example (lAli et al. 2010: B ento et al. 20031: 



Chen et al. 2011: Chen et al. 2011: Chiba et al. 2009: Bento et al. 2006 



Dutta et al. 2011 



Dutta et al. 20ld : 



Dutta et al. 20091 : lHarko et al. 20ld : 



Hirano et al. 2011 



Kuiat et al. 20061 : 1 Novosvadlvi et al. 201Ch 



and references therein) that attempt to explain this accelerated expansion of the universe with stan- 
dard as well as exotic ideas; the canonical and non-canonical scalar fields are the most p romising 
candidates till date. Of late, Robert R. Caldwell and Eric V. Linder (jCaldwell et al. 20051 ) catego- 
rized these scalar field models in two broad classes namely "freezing" and "thawing" dark energy, 
based on asymptotic behavior of the scalar field potential. In thawing models, dark energy equa- 
tion of state wx initially remains at —1 and deviates from —1 near present epoch whereas just the 
opposite behavior of w is witnessed in freezing models. 

In thawing models deviation of wx can lead to w x > — 1 or w x < —1, where w x is value 
of wx at present epoch. For slow-rolling quintessence scalar fields wx moves to w x > —1, 
whereas for slow-rolling phantom scalar fields wx becomes less than —1 at t he present epoch - 
There are other works related to t h ese slow-rolling scalar field thawing mode ls (jDutta et al. 20091 ; 
Gupta et al. 20091 : ISen et al. 201ol : IScherrer et al. 20081 : IScherrer et al. 2008T ) which lead to both 
quintessence and phantom behavior of wx- In these slow-rolling scalar field models with nearly 
flat potential, initially the kinetic energy of the field is much smaller than the potential energy. 
This is because of the initial large Hubble damping which keeps the field nearly frozen at wx = — 1 
at earlier era i.e., in radiation and matter dominated eras. Due to the expansion of the universe, 
energy density of the universe decreases. After the radiation and matter dominated eras, the field 
energy density becomes comparable to the background energy density of the universe resulting in 
the deviation of the field from its frozen state, thereby leading to deviation of wx from —1. 

These slow-roll scalar field thawing models can be characterized by differen t relations be tween 
wx and the scale factor a of the universe. Some typical examples include (ILinder 20081 ) CPL 
parametrization (Eq. (JI])), PNGB models (Eq. ([2])) and algebraic thawing model (Eq. Q). The 
corresponding equation of state parameterizations are given by, 



dwx 
d(ln a) 

dwx 
d(ln a) 



dwx 
d(ln a) 



[l+wx) 3 



: (l + wx) 

F{l + w x ) 
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il) 
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where F is a parameter which is inversely proportional to the symmetry breaking energy scale and 
p and b are two free parameters. 



In this work, we propose a two parameter generalization for this thawing dark energy models 



MS 



dwx 
da 



(l + w x )f(a) 



(4) 



where f(a) is an arbitrary function of scale factor a. In this paper, we have chosen /(a) as 
/(a) = c/a n , where c and n are two otherwise arbitrary param eters. With n = 1 and c = 1 our 
proposal exactly overlaps with CPL thawing dark energy model J Linder 20081 ) . For n = 1 and 1 < 
c < 3, our proposal leads to PNGB thawing dark energy model (ILinder 20081 ). PNGB dark energy 
model is studied ex c lusively for scalar fi e lds dark energy with PNGB potential (IDutta et al. 20071 : 
Frieman et al. 19951 : iKaloper et al. 20061 : iRosenfeld et al. 20071 ) . For suitable choice of the parame- 
ters n and c, our model can approximately reflect algebraic thawing as well. As it turns out, all the 
existing (and probably, upcoming) thawing dark energy models fall in this broad parametrization 
with different values of the parameters n and c. So, rather than proposing individual models, it 
is quite reasonable to construct a generic form of parametrization, analyze it and search for possi- 
ble constraints on the parameters from observations. This is the primary objective of the present 
article. 

We also draw some comparisons among the results obtained for different values of n (i.e., for 
n = 1, n = 1.5 and n = 2) with different fixed values of c and vice-versa. We further provide justi- 
fication for this proposed generalized form of thawing dark energy model against the other existing 
models by comparing them with ours . Moreover, we constra i n our model by latest Supernova Type 



la Data from Union2.1 compilation ( Amanullah et al. 201ol: bavis et al. 2007 : 



Suzuki et al. 2005 



Kowalski et al. 20081: IPerlmutter et al. 19991 : iRiess et al. 19981 : iRiess et al. 2007 



Kessler et al. 2009 



Riess et al. 2009 



Wood-Vasev et al. 20071) , newly released Observational Hubbl e parameter Data 



(jGaztanaga et al. 20091 : Ijimenez et al. 20021 : ISimon et al. 20051: IZhang et al. 20121) . Cosmic Microwave 



Background Shift Parameter Data from WMAP7 year results (jKomatsu et al. 20091 : iKomatsu et al. 2011 



Podariu et al. 200 ll: iRatra et al. 19991) and the latest Baryon Acoustic Oscillati on Data from SDSS 
Data Release 9 ( Percival et al. 201ol : Sanchez et al. 2012 : Samushia et al. 2009 ). For such analyses 
we have five parameters in total namely c, n, w x , and Hq (where $7^ and Hq are matter 
density at present epoch normalized to critical density and Hubble parameter at present epoch 
respectively). Since the value of (the normalized radiation density at present epoch) is very 
low we do not treat it as a parameter and consider = 5.05 x 10 -5 (jBeringer et al. 20120 for 
numericals. Our analysis also helps in comparing the standard diagnostics with model independent 
ones, and reveals the pros and cons of each one. 



The major conclusions of the paper are as follows: 



1. Results obtained for different n values (with different values of c) barely differ from the 
observational point of view. In other words we can say that values of n (with different values 
of c) can hardly affect the matter density parameter O^, ctvera,ge deceleration parameter ^ av 
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and the growth of matter perturbation /. Therefore it is difficult to provide a unique dark 
energy EOS wx for the thawing dark energy models as different values of n with different 
values of c lead to the same cosmological dynamics. 

2. Though the the equation of state parameter does not leave trace on model discrimination for 
thawing dark energy, values of n and c can be constrained by wx — w'x (w' x = ^rjy) plots 



([Caldwell et al. 20051 ) for thawing dark energy models. More importantly, wx — w'x plots 
can also serve as a model discriminator for the thawing dark energy models. The non-linear 
wx — w' x plots can be realized for values of n other than 1 with different values of c. This is 
an important issue as PNGB and CPL parameterizations can result in linear wx — w' x plots. 

3. Most importantly, the model-independent parameters like statefinder pair {r, s} and the so 
called Om3 parameter do play a crucial role in discriminating among different thawing dark 
energy models. Therefore, unlike the standard parameters mentioned in Item 1, these param- 
eters indeed serve as model discriminators for different thawing dark energy models. 



2. The scheme of generalization 

2.1. Generalized thawing dark energy EOS 

We propose a minimal two parameter generalization for thawing dark energy EOS wx as, 

dwx 



da 



(l + w x )f(a) (5) 



where /(a) is an arbitrary function of scale factor a of the universe. We study the dynamical universe 
with radiation, matter and thawing dark energy obeying our proposed EOS wx with /(a) = -4r. 
Our choice of /(a) for generalized thawing model is motivated by the following findings: 

i) for n = 1 and c = 1, our model is exactly same as CPL parametrization (Eq[T]). 

ii) for n = 1 and c = F, our model is exactly same as PNGB model (Eq[2|). 

iii) algebraic thawing case (Eq[3]) can approximately be obtained for certain choice of c and n 
in terms of b and p. 

iv) for any higher values of n, generalized thawing dark energy EOS takes the form 



wx(a) = -1 + (1 + w x ) exp 



c - ,(!-«)> 



(n-1) 

where w x is the value of wx at the present epoch. Expansion of wx{a) about a = 1 gives 



(6) 



wxip) = w x — c(l + w x ){l — a) + -(1 + w x )(c 2 — cn)(l — a) 2 + higher order terms . (7) 
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Fig. 1. — The left fig. is for w x = —0.9 and the right fig. corresponds to w x = —1.1. Black and 
orange plots are for PNGB and algebraic thawing models for w x = —0.9. The dotted, solid and 
dot-dashed lines are for F = 1, 2, 3 in the PNGB thawing case and for algebraic thawing case they 
are for p = b = l, p = b = 2 and p = b = 6. The blue and green curves are for our generalized 
thawing EOS. Green (dotted, solid, dot-dashed) lines are for n = 1.2 (c = 1, 1.2, 1.5). Similarly 
blue (dotted, solid, dot-dashed) lines ar e for n = 1.5 (c = 1, 1.2, 1.5). The the area between solid 



red lines is the allowed thawing region (jCaldwell et al. 20051 ). 



In order to test the validity of our generalized model we show in Fig. [TJ the theoretically 



predicted wx 



X \ W X - dMa) i 

of n and c (we will compare it with direct observational data later in this paper). We find from Fig. 
[H theoretically obtained wx — w' v pla-ne for different combinations of c and n in our model satisfy 
the allowed regions for the same (jCaldwell et al. 20051 ). In Fig. [1] the left plot is for quintessential 
thawing with w x = —0.9 and the right one is for the case of thawing originated in phantom 
scenarios with w x = —1.1. For n = 1 with c = 1 (dotted) we get CPL thawing (Eq. ([1])) and 
for n = 1 with c = 2, 3 (solid and dot-dashed respectively) we get PNGB thawing (Eq. ([2])). The 
plots in black in Fig. [1] indicate these two models in the wx — w' x plane. The orange plots are for 
the algebraic thawing model with p = b = 1 (dotted lines), p = b = 2 (solid lines) and p = b = 6 
(dot-dashed lines). The results with higher values of n are shown by the green (n = 1.2) and blue 
plots (n = 1.5). The dotted, solid and dot-dashed lines in this case corresponds to c = 1, 1.2, 1.5 
respectively. 



dwx 



plots for different thawing models that arise for different values 



2.2. The standard and model independent parameters 



As is well-known, any dark energy model must at least probe three parameters directly from 
observations: 
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i) the present value of equation of state for dark energy (co x ) 

ii) the present value of matter density [Q° M ) 

iii) the Hubble parameter today (Ho). 

Nevertheless, dark energy model building today is tightly constrained by several observations, 
which, taken together, leave out a very narrow window through which the model should pass. So, 
from today's perspectives, apart from the above three good old parameters, the supplementary 
parameters which one needs to address are the following: 



The statefinder {r, s} (jSahni et al. 2003 ) serves as a geometrical diagnostic to probe the prop- 
erties of dark ene rgy in a model independent manner. This pair (r, s} has been studied extensively in 
the earlier works (jTsujikawa 201ol : lAli et al. 201ol : bas et al. 201ll : \u et al. 20091 : IPanotopoulos 20081 ) 
For the late universe (z < 10 4 ), which is well approximated by the presence of matter and dark 
energy, the statefinder pair {r, s} can be expressed as, 

9 

-n x w x (i + w x 



1 + 



1 + -Q, x w x (l + w x ) 



2 H ' 

3 dw x 
-ai ly — - — 

2 da 



1 + w x — - 



1 w x 



1 + w x 



3 w x H 
1 a dw x 



(8) 

(9) 
(10) 

(11) 



3 w x da 

where a is the scale factor of the universe and £l x is the dark energy density parameter and dot 
represents the derivative with respect to cosmic time t. In the late universe we have VL m + £l x = 1, 
O m being the matter density parameter. For ACDM model, it can be checked that the statefinder 
pair {r, s} takes the value r = 1 and s = 0. Any deviation in r from 1 and s from 0, indicates the 
existence of varying dark energy model in the universe. 



The Om parameter proposed by Sahni et al. (jSahni et al. 20081 ). is another tool to distinguish 
the dynamical dark energy from the cosmological constant. The uncertainty in matter density 
parameter allows significant errors in cosmological reconstructions of dark energy. Om parameter 
can in practice differentiate between the models, independe nt of the matter density parameter. The 
Om diagnostic has been studied well in the earlier works (Sahni et al. 20o"3; Nesseris et al. 201C : 
Shafieloo et al. 20121 : Irluang et al. 201 ll : Ishafieloo et al. 201ol : iLu et al. 20081 ). Om parameter is 
defined in terms of Hubble parameter which can directly be measured in cosmological observations. 
The two-point Om diagnostic is given by, 

h 2 (z 2 )-h\ Zl ) 



Om(z 2 ; zi) 



(I + 22) 



(12) 



where h(z) = H(z)/Hq. 



It can be easily seen that for cosmological constant Om(z\ } z 2 ) = a nd when z~\ < zo , 
Om(z\, z 2 ) > (Om(zi, z 2 ) < 0) represents the case of quintessence (phantom) (jSahni et al. 20081 ). 
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This is how Om evaluated at two different redshifts {z\ and z 2 ) can help in distinguishing the 
dark energy model. Needless to mention that this procedure is independent of and Hq. The 
three-point diagnostic Om3 is defined by, 

n i< \ Om(z 2 ;z 1 ) 

Om3{zi,z 2 ,z 3 ) = —— r . (13) 

Om(z 3 ;zi) 

For AC DM model Om3 = 1. 

Another dimensionless parameter, which is useful for determining the beginni ng of cosmic accel 



eration in dark energy model, is the average deceleration parameter q av , defined as (jSahni et al. 20081 ). 

it 2 



1™ = (t 2 - tl ) I q{t)dt ' (M) 



We use Eqs. dU HH [131 USD for evaluating the statefinder pair {r, s}, Om3 and g av for the case 
of our generalization of thawing dark energy model. 

Further more, galaxy cluster formation is not directly influenced by the existence of dark 
energy. The presence of dark energy alters the Hubble expansion rate which affects the growth 
of inhomogeneities in the cold matter sector. In the linear r egime of matter p erturbations, the 
evolution of the inhomogeneities are governed by the relation ( Wang et al. 19981 ) 



d 2 ln5 /d\n5\ 2 lfd\n6\. 1 n M n 3^ 



dihia) 2 \dlna J 2 \dlna J 2 

where the ma tter density contr ast 5 = 5p m /p m with p m the matter density. The growth factor / 
is defined by (IWang et al. 199a ). 

din a 



Eq. (I15h can be written in terms of growth factor / (defined in Eq. (|16h ) as, 



df 
din a 



+ f 2 + -f(i-3w x (i-n m )) 



(17) 



The growth equation can be expressed in terms of the redshift z by the re lation In a = — ln(l + 
z). The growth factor is well approximated by the ansatz (|Wang et al. 1 998) 



/ = n m (z)~< (is) 

where 7 is termed as "growth index". This quantity is affected by dark energy models via £l m (z). 
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3. Compilation of combined datasets 

For the purpose of putting constraints on the generalized thawing dark EOS, we use the 



Davis et al. 2007 


Kessler et al. 2009: Kowalski et al. 2008: 


Perlmutter et al. 


1999: Riess et al. 1998: 


Riess et al. 2007: 


Riess et al. 2009 


Suzu 


a et al. 2005: " 


Wood-Vasev et al. 2007). Observational Hub- 


ble Data (OHD) ( 


Gaztanaea et al. 20091: 


Jimenez et al. 2002: Simon et al. 2005: Zhane et al. 2012). 


Cosmic Microwave Background Data (CMB) from 7 year WMAP results 


(Komatsu et al. 2009: 


Komatsu et al. 2011: Podariu et al. 2001 


: Ratra et al. 1999 


), and Baryon Acoustic Oscillation (BAO) 


Data from SDSS Data Release 9 ( 


Percival et al. 2010: 


Sanchez et al. 2012: 


Samushia et al. 2009). 



There are a total of 607 data points (580 data point from SNe la, 25 from OHD, and 1 each 
from CMB and BAO). We make a combined % 2 analyses of the data sets comprising of all 607 
data points to constrain our model parameters w x , and Hq, as well as to confront with the 
model-independent parameters mentioned in Section 2. This makes our analysis robust. 



3.1. Union 2.1 compilation of Supernova Type la Data 



Luminosity distance {di) measurement of distant supernovae with redshifts z is the first ob- 
servational data to probe the current acceleration of the universe and the dark energy properties 
as well. The most recent compilation of the Supernova Type la Data is given by Union 2.1 dataset. 
The data is tabulated in terms of distance modulus n(z) with redshift z. The distance modulus 
can be written as 

H(z) = 5log 10 (D L (z)) +f i , (19) 

where Dl(z) = HqcIl(z) / c, c being the speed of light and = 42.38 — 51og 10 h with h given by 
H = lOO/iKmSec^Mpc" 1 . 



X 2 of SNe la data is given by, 



E 



Mobs(^) - K z i'i w Xi ^m> H o) 



<Ji 



(20) 



Marginalizing over the nuisance parameter /jlq, one gets the x 2 as : 



j& s (w° x A) = A-B*/C , 



where A, B and C are given by, 

* = E 

i 

B = E 



Mobs(^) - n(znw x , n^no = 0) 



(7; 



(21) 
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C 



(22) 



3.2. Observational Hubble Data (OHD) 

Measurements of Hubble parameters from differential ages of galaxies provide another way 
to probe the late tim e acceleration of the expanding universe. Gemini Deep Deep Survey GDDS 
([Abraham et al. 20041 ). SPICES and VVDS surveys provide the values of the Hubble parame- 
ter at 15 different redshift values and with the dat a release 9 from Sloan Digital Sky Survey 



(SPSS) we have a total of 25 data p oints at present (IGaztanaga et al. 20091 : Ijimenez et al. 20021 : 
Simon et al. 20051 : IZhang et al. 20121 ) . The x 2 function for the analysis of this observational Hubble 



data can be defined as 



15 

E 

8=1 



Hobs(Zi) 



H(zi;w° x ,n° m ,H ) 



(23) 



The latest data release for the variation of Hubble parameter with redshift has been listed in 
Tabled) 



3.3. CMB Shift Parameter Data 



CMB shift parameter R, to a great extent, is a model independent quantity extracted from 
CMB power spectrum. It is given by 



R(z*) = (O m # 2 ) 1/2 / dz/H{i 
J o 



(24) 



where 2* is the redshift value at the time when photons decoupled from matter in the universe, 
can be calculated as (with 0;, being the baryon density parameter) 



z* = 1048[1 + 0.00124(ft fe /i 2 )-°- 738 [l + gi {n m h 2 ) 92 } 
where the functions g\ and <?2 read as 



92 



9l = 0.0783(^/ l 2 )- a238 (l + 39.5(^/i 2 )-°- 763 )- 1 
0.560(1 + 21.1(fi 6 /i 2 ) 1 - 81 )- 1 . 

XcmbKA^o) 



R(z*,w» x ,n» m ,H )-R 



<7R 



(25) 



(26) 
(27) 

(28) 



From WMAP 7 year data (jKomatsu et al. 2009 
we get R = 1.715 ± 0.021. 



Komatsu et al. 201 ll : IPodariu et al. 200J ; iRatra et al. 19991 ) 
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3.4. Baryon Acoustic Oscillation Data 

Baryon Acoustic Oscillation or BAO (observed in the matter power spectrum) is the conse- 
quence of the primordial fluctuations frozen in the large scale structure after the recombination is 
reached. The acoustic peak in the two-point galaxy correlation function was first detected by the 
SDSS by observing the galaxy distribution in the sky and a preferred comoving scale r s , at redshift 
z = 0.35, was measured. 

The dimensionless combination A(z) = Dy (z) ^ Q^H'q / z is independent of Hq and is well 
constrained by the present SDSS da t a. We use the measu r ed value of A(z) at z = 0.57 (A o b s (0.57) = 



0.444 ± 0.014 (jPercival et al. 20ld : ISanchez et al. 2012 ; ISamushia et al. 2003 )) to constrain our 



model parameter space. 

The x 2 f° r the BAO data is defined as 



,o O o x _ [^(0.57) -A(0.57,w° x ,n° m )] 2 



Xt,AO«A) = 0.0162 ' ^ 



3.5. Combined x 2 analyses 

Combining all the datasets from Sections (3.1) - (3.4), comprising of altogether 607 data points, 
the combined x 2 can be evaluated as: 

x 2 (w° x , Q° m , H ) = x| N «, n°J + xhnn«A, H ) + xImb«A, H ) + xIao«A) • 

In what follows, we minimize this combined \ 2 with the observational data sets and search for 
possible consequences by confronting our generalized model directly with observations. 



4. Data analysis and results 

In this section our primary target is to make a combined x 2 analysis for our generalized model 
as proposed in Eq Q with the SNe la, OHD, CMB and BAO data for evaluating the parameters 
space and their la and 2a confidence level (C.L.) limits. We further study these cases to compare 
between the results for n = 1 and other values of n with the different values of c. Our results are 
tabulated in Table [21 El and HI There are five parameters in this generalized thawing model and 
they are n, c, w x , and Ho. We fix the values of n at 1, 1.5, 2 with different values of c so that 
we can compare different thawing models and find the best fit values of other three parameters by 
X 2 analyses. The results of x 2 analyses for PNGB and CPL models are furnished as Case I below 
and the x 2 analyses results for the thawing models with n = 1.5 and n = 2 are given as Case II 
and Case III respectively. 
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4.1. Standard parameters for different values of n and c 

Case I: n = 1 (CPL fc PNGB) 

The x 2 analyses results for n = 1 with different values of c are tabulated in the Table [2j These 
are the cases of PNGB and CPL (c = 1) thawing dark energy models. Here we choose the values 
of c to be 1, 1.5, 2. It is seen from Table that with the increases of the value of parameter c, 
the thawing dark energy EOS transcends from quintessence to phantom domain via cosmological 
constant. During this course of change of EOS (wx), the value of matter density parameter at 
present epoch remains unaltered and present epoch value of the Hubble parameter barely suffers 
any change in its magnitude. Also needless to mention here that the value of total \ 2 remains 
unchanged as is evident from Table [2j 




Fig. 2. — This Fig shows the contours for n = 1. y 2 minimization gives the best fit values which 
are marked as * in the plot. 
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In Fig. [2j the la and 2a contours for n = 1 with different values of c are shown by light 
blue and dark blue shaded regions respectively. The "*" in the plots represents the best fit values 
obtained by x 2 minimization (Table [2]). 

Case II: n = 1.5 

The x 2 minimization results obtained for n = 1.5 with c = 0.5, 1, 1.5 are tabulated in Table[3j 
Here also the best fit results suggest a phase transition from quintessence to phantom nature 
of thawing dark energy EOS as c increases. But the best fit values of present epoch matter 
density parameter remain unchanged. Also one observes that the best fit values of the Hubble 
parameters Hq at the present epoch are hardly changed in these cases. 
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Fig. 3. — This Fig shows the contours for n = 1.5. \ 2 minimization gives the best fit values which 
are marked as * in the plot. 



In Fig. [3l the best fit values (obtained from x 2 minimization) are shown with "*" symbol and 
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the la and 2a contours are given by light blue and dark blue color shadings. 
Case III: n = 2 

The results for n = 2 with c = 0.5, 1, 1.5 are tabulated in Table HI Like the previous two cases 
discussed above, here also one can see the transition from quintessence to phantom behavior of 
thawing dark energy EOS with no change in the present epoch value of matter density parameter 

with the increasing value of c. the present epoch value of the Hubble parameter also barely 
changes. 




Fig. 4. — This Fig shows the contours for n = 2. x 2 minimization gives the best fit values which 
are marked as * in the plot. 

As in the previous two occasions, best fit values and la, 2a contours (obtained from x 2 
minimization) is denoted by "*" and light blue, dark blue color shades respectively in Fig. 01 
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Now we compare different values of n for a particular value of c. For c = 1, one can figure out 
from the Tables [2j [3l [3] that as n value increases from 1 to 2, w x shifts from —0.9992 to —1.0008 
indicating the phase transition from quintessence to phantom nature of thawing. The present value 
of matter density parameter Vfi m remains unchanged and the present value of Hubble parameter Hq 
hardly suffers any change in these cases. The same analogy goes for c = 1.5 case. 

In Fig. the growth of matter perturbation is plotted against redshift z. The upper (lower) 
panel is for f(N = —7) = 0.8 (f(N = 7) = —0.9)). The growth of matter perturbation is identical 
in all the cases suggesting the formation of the same large scale structure in all cases of thawing 
considered here (i.e., for n = 1, n = 1.5 and n = 2). Therefore the growth factor / does not serve 
as a model discriminator but acts as a supplementary probe to confirm correct estimation of cosmic 
structures formed. 
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Fig. 5. — Variation of the growth factor with logarithm of scale factor. The upper (lower) three 
plots are for f(N = —7) = 0.8 (f(N = —7) = 0.9) where N = In a is the number of e-foldings. 
Each plot actually depicts 3 different plots overlapping with each other: for n = 1 with c = 1, 1.5, 2 
and for n = 1.5, n = 2 they are for c = 0.5, 1, 1.5. 

Fig [6] depicts variation of w' x with wx as obtained by diretly using observational data from 
the combined dataset. The plots show that we indeed have non-linear behavior for the generalized 
thawing dark energy model from observations. Comparison of these plots with our theoretical 
predictions, as done in Fig [T] will be interesting. So, Fig [6] goes over Fig [T] which was only a 
theoretical prediction. As it turns out from this figure, the variation of the equation of state indeed 
serves as a model-discriminator for different thawing dark energy models. 
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Fig. 6. — The left plot is for w x > —1 and the right plot is for w x < —1. In both the figs, the blue 
lines are for n = 1.5, the orange line are for n = 2 and the black lines corresponds to n = 1. The 
dotted, solid and dotdashed lines corresponds to c = 1, c = 1.5, c = 2 for the case of n = 1 and 
c = 0.5, c = 1, c = 1.5 for the case of n = 1.5 and n = 2. 
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z 


TT ( \ 

H(z) 






iff — 1 TV If — 1 \ 

(Km. s .Mpc ) 


/ T 7~ 1 TV JT 1 \ 

(Km. s .Mpc J 


0.090 


69 


12 


0.170 


83 


8 


0.270 


77 


14 


0.400 


95 


17 


0.900 


117 


23 


1.300 


168 


17 


1.430 


177 


18 


1.530 


140 


14 


1.750 


202 


40 


0.480 


97 


62 


0.880 


90 


40 


0.179 


75 


4 


0.199 


75 


5 


0.352 


83 


14 


0.593 


104 


13 


0.680 


92 


8 


0.781 


105 


12 


0.875 


125 


17 


1.037 


154 


20 


0.24 


79.69 


3.32 


0.43 


86.45 


3.27 


0.07 


69.0 


19.6 


0.12 


68.6 


26.2 


0.20 


72.9 


29.6 


0.28 


88.8 


36.6 



Table 1: Hubble parameter (H(z)) versus redshift (z) data from (jGaztanaga et al. 20091 ; 
Jimenez et al. 20021 ; ISimon et al. 20051 ; IZhang et al. 20121 ) . Where H(z) and an are in km ■sT 1 
MpcT 1 . 
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n 


c 


best-fit values of 

(w , n^, h ) 


Minimum 
value of x 2 


1 


1 


(-0.9992, 0.279, 70.47) 


576.26 


1 


1.5 


(-1.0008, 0.279, 70.49) 


576.26 


1 


2 


(-1.0023, 0.279, 70.50) 


576.26 



Table 2: Best-fit values of parameters and minimum values of x 2 from combined x 2 analyses of 
SNIa, CMB, OHD and BAO data for Case I. 



n 


c 


best-fit values of 

(w , n^, h ) 


Minimum 
value of x 2 


1.5 


0.5 


(-0.9979, 0.279, 70.46) 


576.26 


1.5 


1 


(-0.9999, 0.279, 70.48) 


576.26 


1.5 


1.5 


(-1.0017, 0.279, 70.49) 


576.26 



Table 3: Best-fit values of parameters and minimum values of x 2 from combined x 2 analyses of 
SNIa, CMB, OHD and BAO data for Case II. 



n 


c 


best-fit values of 
{wq, n° m , H ) 


Minimum 
value of x 2 


2 


0.5 


(-0.9985, 0.279, 70.47) 


576.26 


2 


1 


(-1.0008, 0.279, 70.49) 


576.26 


2 


1.5 


(-1.0027, 0.279, 70.50) 


576.26 



Table 4: Best-fit values of parameters and minimum values of x 2 from combined x 2 analyses of 
SNIa, CMB and OHD and BAO data for Case III. 
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4.2. Model-independent diagnostics 



In Fig. [7] we show the variation of the statefinder parameters {r, s} with redshift z for n = 1 
case which is known as CPL for c = 1 or PNGB for other values of c. The dashed, dotdashed and 
dotted plots are for c = 1, 1.5, 2 respectively. This plots bear the clear signature of thawing as one 
can see that for higher values of z, the statefinder r tends to 1 and the statefinder s to 0. This 
is because wx = —I as z increases and since in present epoch wx deviates from — 1, r and s also 
deviates from 1 and respectively. The same features are also observed in the cases of n = 1.5 and 
n = 2 in Fig. [8] and Fig. [9] respectively. 



Fig. 7. — Statefinder r and s as a function of redshift z for n=l and c=l (dashed), 1.5 (dotdashed), 
2 (dotted). 




0.0 0.5 1.0 




Fig. 8. — Statefinder r and s as a function of redshift z for n=1.5 and c=0.5 (dashed), 1 (dotdashed), 
1.5 (dotted). 



In Fig. [TUl we show the variation of the Om3 parameters with the redshift Z3 while Z\ and zi 
are kept at Z\ = 0.2 and Z2 = 0.57. The plot at the extreme left of Fig. [10] shows this variation 
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0.0 0.5 1.0 



0.0 0.5 1.0 1.5 2.0 



Fig. 9. — Statefinder r and s as a function of redshift z for n=2 and c=0.5 (dashed), 1 (dotdashed), 
1.5 (dotted). 



with n = 1 for c = 1 (dashed), c = 1.5 (dotdashed) and c = 2 (dotted). The middle and the right 
plots of Fig. [10] show similar variations for n = 1.5 and n = 2 respectively with c = 0.5 (dashed), 
c = 1 (dotdashed) and c = 1.5 (dotted). As Om3 is a three point diagnostic, we need three redshif t 
points to meas ure its value. We fix two redshift points z\ and Z2 with z\ = 0.2 ([Blake et al. 201 ll ) 
and Z2 = 0.57 ([Sanchez et al. 2012 ) and allow z% to be a variable. All the variation starts from a 
point where z% = Z2 that leaves Om3 = 1 and the immediate deviation of Om3 from 1 suggests 
the varying nature of dark energy which is indeed the nature of variation in thawing dark energy 
scenario in this case. 
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Fig. 10. — The left plot shows the variation of Om3 parameter as a function af redshift z^ with 
z\ = 0.2 and zi = 0.57 for n = 1 with c = 1 (dashed), c = 1.5 (dotdashed), c = 2 (dotted). The 
middle and the right plots shows the same for n = 1.5 and n = 2 respectively with c = 0.5 (dashed), 
c = 1 (dotdashed), c = 1.5 (dotted). 



In Fig. [TTJ the average deceleration parameter q av has been plotted. It is seen from the plots 
that all of them overlap with each other. It is thus evident that average deceleration parameter 
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is not capable of being a model discriminator, but it does indicate the transition period from the 
deceleration to acceleration phase. In this case this transition occurs between the redshift range 
6.25 < z < 6.5 as is clear from the best fit plots. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 



Fig. 11. — q av vs z plot for c=l and n=l, 1.5, 2. 



5. Discussions and Conclusions 

In this work, we proposed a two parameter generalized EOS, wx, for thawing dark energy 
models and studied the dynamics of spatially flat FRW universe containing radiation, matter and 
dynamical dark energy. This proposal of ours is a minimal generalization of thawing dark energy 
EOS and is given by, 

w ' x (a) = (l + w x )^ . (30) 
ar 

This leads to wxifi) = —1 + (1 + w x )a c for n = 1 and for other values of n, wx{a) = —1 + (1 + 
w x ) exp[ ^ n ^ (1 — a( 1-n ))], where each value of the parameters n and c defines a thawing model. 
Here w x is the value of wx at the present epoch. As we have demonstrated, our model goes over the 
well-known thawing dark energy models such as CPL, PNGB and algebraic thawing, for suitable 
choice of the parameters n and c. 

We have elaborately discussed the cases with n=l, n=1.5 and n=2 for different values of c 
(c = 1,1.5,2 for n = 1 and c = 0.5,1,1.5 for n = 1.5, 2). We have shown that though the 
parameter c is very important for the slope of the wx{o) vs. a plot, it barely changes the dynamics 
of the universe. This is quite evident from the average deceleration parameter (g av (z)) vs redshift 
z plot (Fig. [TTj) . growth parameter plots (Fig. [5]) etc and also from the present values of matter 
density parameter fiJri and the Hubble parameter Hq as well. In this context it is therefore very 
important to mention that fine tuning of c does not, at all, effectively change the observables like 
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the values of density parameters and the Hubble parameter at the present epoch (which barely 
suffers any change as we move from one model to another by changing the n and c values (vide 
Tables [2j [3] and HJ). except for the variation of wx{o) with the evaluation of scale factor a of 
the universe. Here it is necessary to mention that in spite of treating the present epoch value of 
radiation density param eter 0° as a paramet er in the numerical analysis, we have chosen its value 



to be = 5.05 x 10" 5 (|Beringer et al. 2012 1. This is because of the small value of which will 
not change the total density parameter upto four decimal places and therefore considering it as a 
parameter will not affect density parameters $7^ or $7^ upto four decimal places. 

Also we would like to conclude that different values of n would lead to same cosmological 
dynamics for a particular value of c which is evident from average deceleration parameter plot (Fig. 
[TTj) and growth of matter perturbation plots (Fig. [5]). These plots clearly demonstrate that it is 
hardly possible to distinguish between the results the results for different thawing models (related 
to different values of n and c) . But the statefinder pair r, s and Om3 parameter can discriminate 
between the thawing models. These are shown in Figs. El El EJ [TTJl 

In this context we have shown in Fig. [6] that wx — w' x plots can, as well, serve as another 
discriminator of these thawing models. Nevertheless, it is shown in Fig. [6] that wx — w' x plots are 
non-linear for n = 1.5 and n = 2. For the existing thawing models (e.g., PN GB and CPL cases), 



wx — w' x plots are strictly linear. In a recent work, Ali et al. (|Ali et al. 20091 ) have found this kind 
of nonlinear w — w' plots arising from scalar field models. So our generalization can also produce 
them for values of n other than 1. Moreover from Fig. [6l it can be noted that as n takes higher 
values, only lower values of c are allowed for thawing dark energy models. 

As mentioned, this is a minimal generalization with the two parameters c and n and one 
boundary condition given by wx(z = 0) = w x , z being the redshift. There may exist other 
generalization with more than two parameters. So selection can be made on the basis of Akaike 
Information Criterion (AIC) and the Bayesian Information Criterion (BIC) that are defined as, 

AIC = -21n(£) + 2p, (31) 
BIC = -21n(£)+plniV , (32) 

where C is the maximum likelihood value which is given by exp(— x^ in /2), p is the number of 
model parameters and N is the number of data points used to find the minimum value of the 
X 2 denoted by Xmin- We show the A AIC and ABIC values in Table [5j Usually from statistical 
analysis, it is inferred that the models having ABIC in the range — 2 are strongly supported, 
models with ABIC > 2 are moderately supported, and those with ABIC > 6 are unsupported 
from perspective of a given data. However, in cosmology, with the rapid increase of the number 
of data points (we remind the reader that we have used combined dataset), Eq ([32]) shows 
that the ABIC value is always going to increase with addition of a model parameter p. This 
does not essentially mean that the model with least number of parameters are always favored by 
observations, though it may appear to be so. For example, we know ACDM model (with the 



least number of parameters) ()Li et al. 20101 ) fits the SNe la data only in the low redshift region i.e, 
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Model 


AAIC 


ABIC 


CPL 








PNGB 


2 


6.39 


Our Model 


4 


12.78 



Table 5: A comparative analysis of the values of Information criteria using combined x 2 analyses 
of SNIa, CMB, OHD and BAO data. 



for z « 1, and in this vein, most of the models pay the price just because they have additional 
parameters, though they, i n fact fair well with o bservations. On this note it should be mentioned 
that, as demonstrated in ([Parkinson et al. 20061 ) the above information criteria should better be 
replaced by Bayesian Evidence calculation, which gives a value after integrating over all probable 
states, and hence, does not suffer from any such limitations of AIC or BIC. Hence, nowadays, most 
of the cosmological models are relying more on Bayesian Evidence calculation, rather than AAIC 
or ABIC calculation. We hope to address this issue in near future. 
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